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^^Abstract 

We derive the exact time-evolution for a general quantum system under the influence of pure phase noise and demonstrate that 
for a Gaussian initial state of the bath, the exact result can be obtained also within a perturbative time-local master equation 
approach already in second order of the system-bath coupling strength. We reveal that this equivalence holds if the initial state 
' _ . '6f the bath can be mapped to a Gaussian phase-space distribution function. Moreover, we discuss the relation to the standard 
j3Ioch-Redfield approach. 
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1. Introduction 

The coherent evolution of small quantum systems is typ- 
ically influenced by its interaction with environmental de- 
grees of freedom, which results in quantum dissipation and 
decoherence. These ubiquitous phenomena play a crucial 
role in various fields of physics and chemistry, such as quan- 
tum optics, electron transfer reactions [1,2], the electron 
transport through molecular wires [3,4], and in particu- 
lar in quantum information processing, where we recently 
witnessed significant experimental progress [5-7]. The op- 
timisation of the coherence properties of quantum devices 
certainly requires a good theoretical understanding of the 
processes that induce decoherence. 

The environment of a quantum system is frequently mod- 
elled as an ensemble of harmonic oscillators that couple to 
the system [8-11]. If the coupling is linear in the oscillator 
position, one can formally eliminate the environment to ob- 
tain a closed equation for the dissipative quantum system. 
Such equations are in general not easy to deal with and, 
accordingly, only a few exact solutions exist in dissipative 
quantum mechanics, e.g. for the dissipative harmonic os- 
cillator [12-15] and its parametrically driven version [16]. 
Recently, an exact solution has been found also for the dis- 
sipative Landau-Zener problem at zero temperature [17]. A 
whole class of system-bath models that can be solved ex- 
actly are those in which the system Hamiltonian and the 
system-bath coupling commute [18-27]. Herein we focus on 



such so-called pure phase-noise models. 

Even though pure phase noise allows an exact solution 
of the reduced quantum dynamics, it is sometimes conve- 
nient to employ a perturbative master equation approach, 
such as the Bloch-Redfield equation [28] . On the one hand, 
those approaches may provide direct access to the dephas- 
ing rates avoiding tedious algebra and in particular in the 
limit of weak system-bath coupling, they are expected to 
give quantitatively good results. On the other hand, it is 
possible to test their quality when an exact solution is avail- 
able. We here show that for phase noise, the results of ap- 
proximate master equations can even be exact, despite the 
fact that they are based on second-order perturbation the- 
ory 

After introducing our system-bath model in Section 2, 
we present in Sections 3 and 4 a master equation for weak 
system-bath coupling and the exact time-evolution of the 
reduced system density operator, respectively. In Section 5, 
we argue why the exact solution complies with the master 
equation. 



2. System-bath model 

We model the dissipative quantum system by coupling 
the central system to a quantum heat bath that consists of 
harmonic oscillators, so that the total system-bath Hamil- 
tonian reads [10, 11] 
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H = H, + H,^+H^, (1) 

and Ho = Hs + describes the system and the bath in 
the absence of the coupHng. While we will not specify the 
system Hamiltonian, we employ for the bath an ensemble 
of independent harmonic oscillators with the Hamiltonian 



(2) 



Here, 5^, and bk are the usual creation and annihilation 
operators which obey the commutation relation [bk, = 
^fefc' • We assume that the system couples linearly to the bath 
with a hcrmitian system operator X, so that the interaction 
Hamiltonian reads 

ifsb = nx^, (3) 

with the effective bath coordinate 



(4) 



If the coupling operator X commutes with the system 
Hamiltonian, [X, Hg] = 0, the coupling (3) does not induce 

transitions between system cigcnstatcs and, thus, consti- 
tutes pure phase noise. Henceforth we shall focus on this 
case. 

We choose an initial condition of the Feynman- Vernon 
type, i.e. one for which the total density operator R at 
time t = can be factorised into a system and a bath 
contribution p and pb, respectively, i.e. 



i?(0)=p(0)pb(0). 



(5) 



The bath itself is frequently assumed to be initially at ther- 
mal equilibrium. However, if the initial expectation value of 
the coupling operator X does not vanish, the coupling (3) 
entails a force on the bath oscillators. Then the natural ini- 
tial state pb = Pb(0) of the bath is rather a displaced ther- 
mal state which falls in the class of non-squeezed Gaussian 
states. A convenient basis for these states is provided by 
the coherent states {\(3k)} defined by the eigenvalue equa- 
tion 6^,1/3^,) = /3fc|/3fe). Owing to the overcompleteness of this 
basis, any hermitian operator can be written in a diagonal 
form, which assigns to each operator a P-function [29,30]. 
In particular, the bath density operator can be written as 



(6) 



where <P(3k denotes integration over the complex plane. 
Henceforth, we assume that the P-function of each oscilla- 
tor fc is a Gauss function, such that 



Pk{f3k, Pk) = exp ^ '-^ '— 

wnk \ rik 



(7) 



As a central property of a Gaussian state, all expecta- 
tion values are fully determined by Pk = {bk)h and Uk = 
{bkbk)h - \{bk)b\'^, where 

(...)b=trb(pb...) (8) 

denotes the expectation value with respect to the bath 
state pb- An important particular case is the canonical 



ensemble of the bath at temperature T such that pb oc 
exp(—iJb/fcBr), which corresponds to {bk)h = and 2nk = 
coi\i{tiWk/2k^T) - 1. 

The dynamics of the system plus the bath is governed by 
the Liouville-von Neumann equation 



R{t) = CR{t) , 



(9) 



with the Liouvillian C{. . .) = —i[Hsh{t), . . .j/h. The tilde 
denotes the interaction-picture representation with respect 
to Hq, i.e. A{t) = UQ{t)AUo{t), with the free propagator 
Uo{t) = exp(— iTJoi/fi-)- The interaction-picture represen- 
tation of the effective bath coordinate ^ is easily obtained 
from bk{t) = 6fe exp(— iojfct), while for pure phase noise, 
X{t) ^ X, owing to the relation [X, Hs] = 0. For the same 
reason, Hs and X possess a complete set of common eigen- 
states for which the respective eigenvalues are de- 

noted by En and X^. 

We are exclusively interested in the state of the system, 
so our goal is to find the reduced density operator p{t) = 
tibRit). In the subsequent sections, we derive explicit ex- 
pressions for the reduced dynamics. 

3. Master equation approach 

A common and successful approach to dissipative quan- 
tum dynamics is provided by master equations, i.e. dif- 
ferential or integro-differential equations of motion for 
the reduced density operator [30-32]. There exist vari- 
ous formally exact quantum master equations in time- 
convolutionless [16, 33^2] and time-nonlocal form [43-45] 
which, however, generally cannot be solved explicitly and, 
thus, one often has to resort to a perturbativc treatment. 
The need for approximations can stem from the fact that 
the model itself cannot be solved exactly, or from a tech- 
nically difficult structure of the master equation, or even 
both. 

Here, we employ a time-convolutionless quantum master 
equation of the form 



m = mm ■. 



(10) 



with a time-dependent superoperator K,{t). Note that there 

arises no inhomogcneity since wc arc assuming a factorising 
initial state of the form (5), which leads to a linear equa- 
tion of motion [46]. Equation (10) is formally exact and 
possesses an apparently simple form, but it generally can- 
not be solved analytically. Thus, it is convenient to expand 
the generator K,{t) in powers of the interaction, i.e. K{t) = 
ICf {t). By a direct calculation [47,48] or by using a time 
convolutionless projection operator technique [39,49,50] it 
is possible to obtain an expression for the ^th order gen- 
erator JCf . In doing so, we assign to a superoperator Q 
of the total system a reduced superoperator {Q) defined 
by its action on a system operator Y, that is {Q)Y = 
trb{^(^Pb)}- With this notation, the time-convolutionless 
generators read ICi{t) = {C{t)), and 
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ICeit) 



dti / dt2 
Jo 



(11) 



for £ = 2, 3, . . . The symbol ((. . .})oc denotes an ordered 
cumulant [31,48], i.e. a sum of certain products of reduced 
superoperators of the form {C{t)C{ti) . . .). 

The fact that the Lionvilhans at different times generally 
do not commute makes it practically impossible to write 
down an explicit expression for the iih cumulant for large L 
However, for weak system-bath coupling it is possible to ne- 
glect higher than second order terms in the expansion of the 
generator, i.e. we may approximate lC{t) « lCi{t) + A^2(0- 
Fortunately the second time-ordered cumulant takes the 
simple form {{C{i) C{h))) = {C{t)C{h)) - {C{i)) {C{h)) . 
Considering now explicitly the interaction Hamiltonian (3), 
we obtain the "standard" time-local weak-coupling equa- 
tion 
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m= -mt)uxrp{t)\ 

' dT(s(t,t-T)[x,[xrp{m (12) 

+\A{t,t-T)[x,{x,pmi), 

where {A, B} = AB + BA denotes the anti-commutator, 

and where we have defined the symmetric and anti- 
symmetric correlation functions 



5(t,i') = ^({A|(t),A|(i')})b, 
A{t,t') = \m{t),m')]%, 



(13) 
(14) 



of the operator- valued fluctuation A^(f) = ^(i) — (C(O)b- 
Note that for a Gaussian initial state pb of the bath, the 
mean value (C(O)b of the bath coordinate in general does 
not vanish. Thus, it explicitly appears in the master equa- 
tion (12). 

It is convenient to expand the master equation (12) into 
the eigenbasis of the system-bath interaction. We then ob- 
tain for a matrix element {m\p\n) = pmn the differential 
equation 

Pmn{t) = [ - l{Xm - ^n)(e>))b 

- (X„ -Xnf dr Sit, t - t) (15) 

- i{Xl - Xl) [ dr Ait, t - t)1 pmnit) . 

Jo J 

For the diagonal matrix elements pnn the right-hand side 

of Eq. (15) vanishes, i.e. the populations remain constant 
in time as one expects for a pure dephasing model. This 
implies that generally the system will not reach thermal 
equilibrium. Although no energy is exchanged and, thus, 
the interaction is dissipationless, the relative phases be- 
tween eigenstates will be randomised. As a consequence, 
off-diagonal elements of the reduced density matrix may 
decay, which reflects the process of decoherence. 



Using the coherent state representation (7) for a Gaus- 
sian bath state, we can evaluate the mean value of the bath 
coordinate ^ and the correlation functions in an explicit 
form and find 



m)h = 2^ (gA 6-"^'=* +glPt 6-'=*) (16) 

fe 

Sit,t - t) = Y,\9k\'' cosicOkT)il + 2nk) (17) 
fc 

Ait,t-T) = -Y,\9kfsmi^kT). (18) 

k 

Note that the correlation functions Sit, t') and Ait, t') de- 
pend on time differences t — t' , only, with Ait, t') inde- 
pendent of the initial state of the bath. If the correlation 
functions vanish sufficiently fast, it is possible to employ a 
Markov approximation, i.e. to extend the upper integration 
limits in Eq. (15) to infinity. Then, the master equation (15) 
reduces to the standard Bloch-Redfield equation [28] and, 
moreover, is of Lindblad form, so that the complete posi- 
tivity of the reduced density operator is conserved even for 
arbitrarily short times [51,52]. 

4. Exact solution 

The dynamics of a system subject to pure phase noise 
can, at least in principle, be solved analytically [18-27] . The 

formal solution of the Liouville-von Neumann equation (9) 
after tracing over the bath's degrees of freedom reads 

pit)=XY^{Uit)RiO)U\t)}, (19) 

with the propagator 



Uit)=T exp 



drH^^ir) 



(20) 



and the time-ordering operator T. Although the considera- 
tion of the time ordering can often be quite cumbersome, it 
nevertheless can be accomplished for the model discussed 
here. We deferred the explicit derivation to Appendix A, 
where we find for the elements of the reduced density ma- 
trix the exact expression 

Pmnit) = p„„(0)e'(^-^")^«-(^'"-^")^E. 

X n Xk{zkit)[Xm - Xn], zmiXm - X„]} , 



with the time-dependent phase 

\9k[' 



(21) 



(22) 



(pit) = — ^ [ukt - siniukt)] , 
k '^fe 

the quantum characteristic function [30] 

XkiXk,K) = le^'''''^-^'^0^Pkif3k,f3l)d^f3k, (23) 

and the time-dependent complex number z^it) = £(^[1 — 
expiiu}kt)]/cOk- To arrive at Eq. (21), we only used the as- 
sumption that initially the system and the bath are uncor- 
related and that the initial bath state ph factorises with 



3 



respect to the modes k [sec Eqs. (5) and (6)]. For the case 
of a Gaussian distribution of the bath modes, it is possible 
to calculate the integral in the characteristic function (23) 
explicitly. We finally find the exact time evolution of the 
reduced matrix element 



(24) 



Pmn{t) = Pmn(O) exp|-(X„ - XnfA{t) 

+ i[{Xl - X^)<f>{t) + (Xm - Xr^Mt)] } , 

with the phases (j){t) defined in Eq. (22) and 

^(t) =2^Im(^[l-e'^'=*]) . (25) 



The time dependent damping amplitude A{t) does not de- 
pend on the mean values pk of the bath modes and reads 



k 



\9k\ 



. 1 — COs(Wfct) 



(1 + 2nk) . (26) 



Note that a similar result was obtained recently for a bath 
initially in a squeezed thermal state [53]. 

Upon computing the time-derivative of the exact solu- 
tion (24) and noting that the relations 



f 
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dTA{t,t-T) 



dTS{t,t- t) 



(27) 
(28) 

(29) 



hold, we find the surprising fact that the exact solution 
obeys the quantum master equation (15)! Or put differ- 
ently, for pure phase noise, the exact result can be obtained 
within second-order perturbation theory from the master 
equation (12). For large times t, this master equation be- 
comes the standard Markovian Bloch-Redfield equation. 
Thus, we find that the latter contains the exact long-time 
limit of the rates (27)-(29). 

Before disciissing the relation between both approaches 
in more detail, we like to close this section by writing the 
exact solution (24) also in terms of the usual bath spectral 
density [10,11] 



J{uj) = ^ \gkfS{u) - oJk) ■ 



(30) 



For the important special case of a heat bath that is initially 
at thermal equilibrium, we find (fi{t) = 0, while the phase 
(j){t) and the damping amplitude A{t) read 



(j){t) = y dw J{u> 



A(i) = j dw J(w)- 



ut — sm{wt) 
■ cos{u!t) 



coth I 



\2kT,TJ 



(31) 
(32) 



5. When second order is exact 



We have seen that the time-local master equation (10) 
derived within second order perturbation theory provides 



the exact time evolution of the reduced density matrix, 
which implies that in the expansion of the Liouvillian K,{t), 
all higher order contributions vanish. This on the one hand 
nicely simplifies practical calculations. On the other hand, 
it poses the question whether we face a coincidence or 
whether there is any profound reason for the equivalence. In 
order to underline the latter point of view, we now demon- 
strate that for phase noise, the time-ordered cuniulant in 
the £th order generator (11) is proportional to the usual 
classical cumulant of the initial bath state. Consequently, 
wc! can argiie that for the Gaussian initial state (7), the 
series lC{t) = J2e^ii^) terminates after £ = 2, which im- 
plies that the second-order time-local master equation (12) 
is exact. 

We start out by defining averages with respect to the P- 
function as 



J ■■■l[P{Pk,P*k)d^fik. (33) 



With this notation, the exact solution (21) reads 



Pmn(t) = ((ixpl / drfmnir 



Pmn(O), (34) 



with the complex valued function 

fmnit) = i{Xl - Xl)cj>{t) - ^ {(X„ - X„)^ 

fe 

-{X^-Xn)[zk{t)Pl-zl{t)pk]]. 

(35) 

The average in Eq. (34) is obtained from a distribution 

function for the c- numbers /3fe. Thus, it can be formally 
considered as the averaged solution of a stochastic differen- 
tial equation that obeys a time-local differential equation 
of the form (10), but with the generator K, now being a c- 
number, not an operator. Thus, we can adapt the line of 
argumentation given by van Kampen for classical Gaussian 
stochastic processes [48] : Differentiating the Taylor expan- 
sion of Eq. (34), we find 

{fmn{t))p+ I dti{fmn{t)fmn{tl))p 

Jo 



dti / dt2 {fmn{t)fmn{tl)fmn{t2)) P + ■ ■■ 

Jo 



n(0). 

(36) 

A time- local equation of motion for Pmn(t) can be obtained 
by inserting pmn{Q) from Eq. (34), which yields [48] 

Pmn(t)= dt2... dte-l 

X {{fmn{t)fmn{tl) ■ ■ ■ fmn{tl-l)))p Pmn{t) , 

(37) 

where ((. . .))p denotes the cumulants with respect to the P- 
function. Note that for the cumulants of a classical process, 

time-ordering is not relevant [48] . Thus the only difference 
of this expansion and the one in Eq. (11) for the quantum 



4 



master equation is that the latter contains time-ordered 
cumulants. 

For a Gaussian P- function, all cumulants of and jS^ be- 
yond second order vanish [54] . Since fmn (t) is linear in these 
variables, the same is true for the cumulants in Eq. (37) 
and, consequently, only the terms with £=1,2 contribute 
to this expansion. Evaluating the expansion coefficients ex- 
plicitly, one finds that they arc identical to those of the 
second-order time-local master equation (15). 

The equivalence of the second-order master equation and 
the exact solution is based on two requirements: First, the 
coupling operator X needs to be diagonal in the eigenba- 
sis of the system, so that its^interaction-picture represen- 
tation is time-independent, X(t) — X, and, thus, it can be 
effectively treated as a c-number. Hence the quantum me- 
chanical time-ordering affects only the bath coordinate ^ (t) 
for which wc can express multi-time expectation values as 
cumulants of the P-function. In that way, we can circum- 
vent the tedious task of normal-ordering the operators bk 
and bl.. With this precondition, secondly, the Gaussian ini- 
tial state of the bath ensures that the cumulant expansion 
terminates after the second order and agrees with the ex- 
pansion of the master equation (10). For any non-Gaussian 
state, infinitely many higher-order cumulants are non-zero 
both in the classical case [55-57] and in the quantum me- 
chanical case [58-60]. Consequently, the expansion of the 
Liouvillian is of infinite order and any truncation represents 
an approximation. 

Let us finally stress that the second-order Nakajima- 
Zwanzig master equation [43,44], which was not considered 
in this work, can be expressed in terms of cumulants (the 
so-called partial cumulants), as well [61,62]. Note that the 
ordered and the partial cumulants up to second order coin- 
cide. However, the second-order Nakajima-Zwanzig equa- 
tion is not of a time local form and therefore cannot yield 
the exact result for the model discussed here [63]. Thus, 
phase noise constitutes an example for which the time- 
local approach is more accurate than the time-nonlocal one 
when comparing their perturbation expansions up to the 
same order. This outcome is in agreement with some re- 
cent findings for harmonic oscillator baths [64, 65] and for 
spin baths [66,67]. For example for a two- level system cou- 
pled via XY Heisenberg interaction to a spin bath, the dif- 
ferences of both approaches have been analysed quantita- 
tively [66]. Nevertheless, we do not give a general recom- 
mendation in favour of one or the other approach because 
the quality of each seems to be model dependent [6 1 , 68-70] . 

6. Conclusion 

Quantum systems under the influence of pure phase noise 

represent an important special case of dissipative quan- 
tum mechanics owing to the existence of an exact solution. 
Moreover, on short time scales, on which the coherent sys- 
tem dynamics cannot manifest itself, the behaviour of the 
phase noise model is even generic [23]. Here, we have pre- 



sented the explicit exact solution for a quantum system 
under the influence of phase noise with a general Gaus- 
sian initial bath state. Thereby, we have demonstrated that 
the coherence decay is determined by the symmetric bath 
correlation function, while the anti-symmetric correlation 
function gives rise to a time-dependent phase shift. In turn, 
from the exact relations (27)-(29), one can obtain informa- 
tion about the spectral properties of the bath by comparing 
our results with the experimentally observed dephasing at 
short times. 

Despite the exact solvability of the phase noise prob- 
lem, it is often convenient to study the resulting dephasing 
within a master equation approach based on second-order 
perturbation theory in the system-bath coupling. For the 
(time-nonlocal) Nakajima-Zwanzig equation, this consti- 
tutes an approximation. For the time-local version of such 
a master equation, by contrast, we have found that it pro- 
vides the exact solution. After noticing that this facilitates 
practical calculations, one might wonder why and when this 
equivalence holds true. By mapping the initial bath den- 
sity operator to a P-function, we showed that a formal ex- 
pansion of the time-local master equation for phase noise 
is in fact an ordinary cumulant expansion. Consequently, 
for a Gaussian initial bath state, all terms beyond the sec- 
ond order vanish, so that the master equation becomes ex- 
act. Thus for a bosonic heat bath, there are two conditions 
for the exact agreement: First, the system-bath coupling 
must commute with the system Hamiltonian constituting 
the case of pure phase noise and, second, the initial state 
of the bath must correspond to a Gaussian P-function. If 
one of these conditions is violated, there might still exist 
an exact solution, but it can no longer be obtained within 
second-order perturbation theory, as for example is the case 
for the dissipative harmonic oscillator [13, 14]. 

The second-order time-local master equation employed 
in this work agrees with the exact solution at any time. In 
particular in the long-time limit, it becomes Markovian and 
identical to the standard Bloch-Redfield master equation, 
which for pure phase noise is of Lindblad form. This also 
explains the previously observed "excellent agreement" [71] 
between the exact dynamics and results obtained within 
Bloch-Redfield theory. 
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Appendix A. Derivation of the exact solution 

In this appendix, we outline the derivation of the exact 
reduced dynamics discussed in Sec. 4. As a first step, we 

perform a transformation to the interaction picture with re- 
spect to Hq, so that the coupling Hamiltonian (3) becomes 
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H,k{t) = Vit) + vHt), 



(A.l) 



with V{t) = gfc6fc exp(— iwfet) and the hermitian 

system operator X. The two operators V and do not 
commute. However, their commutator is an operator in 
the Hilbert space of the system, while being a c-number 
in the bath Hilbert space. To be specific, one obtains 
[V{t),V^{t')] = fit- t'), where 



(A.2) 



Hence, we can use the Baker-Campbell-Hausdorff formula 

(see Rcf . [30] ) to express the time-ordered propagator (20) 
as a product of two exponentials 

C/(i)=exp|l^*drH,b(T)| 

xexp|-^^*dr^"dr7(T-T') (A.3) 

x[e{T-T')-e{r'-T)]y 
The first exponential in Eq. (A.3) can be written as 

expjl ^ dri?,b(r)| = ^/^^(zfcX) (A.4) 

with the complex number Zk = ^^[1 — exp{iukt)]/uk and 
the displacement operators Dk{Y) = exp(F&| — V^bk). 
The second exponential in Eq. (A.3) provides the time- 
dependent phase factor ex.p{iX^ (p{t)} with 

<^(t)= / dr / dr' V|<?fepsin(a;fer'). (A.5) 
Jo Jo J, 

We will now consider the elements of the reduced density 
matrix in the eigenbasis {|n)} of the system Hamiltonian. 
Using the initial condition (5), Eq. (19) becomes 

Pmn{t)= trb(m|[/(t)p(0)pbC/t(t)|n) (A.6) 
= P™n(0)e'[^^-^n]^W 

l[Dl{X^Zk)l[Dl,{X^Zk')) , (A.7) 



where (. . .)b = trb(/fb---)- -^Y virtue of the relations 

DliY) = Dli-Y) and 

Dk{Y)Dk{Z) - Dk{Y + Z) cxp{{YZ^ - Y^Z)/2} , (A.8) 

which hold for any commuting system operators Y and Z, 
we obtain 



'l[Dkizk[X^-Xn]) 



(A.9) 



Note that an additional phase rjmn = 2 X)fe kfe |^ ^^{^m-^n) 
vanishes, since the system operator X is hermitian. 

It remains to evaluate the expectation value in the sec- 
ond line of Eq. (A.9). This is readily established by writ- 
ing the bath state pb in its P-function representation [see 



Eq. (6)] and noticing that expectation values of normal or- 
dered products of annihilation and creation operators are 
identical to the moments of the P- function, where and bk 
have to be replaced by PI and Pk, respectively [30]. Thus, 
by using the Baker-Campbell-Hausdorff formula 

exp(A6^ - Xlbk) = exp(A6^) exp(A*6fe) exp(-|A| V2) 

(A.IO) 

for each mode k, we write the second line of Eq. (A.9) in 
its normal-ordered form 

(A.ll) 

where {■ . ■)p denotes the average (33) with respect to the 
P-function. Using also the fact that the P-function (7) fac- 
torises with respect to the modes k, we arrive at Eq. (21). 
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